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This paper proposes an accurate model of Photovoltaic Module. The model is 
presented in details and simulated by Matlab Simulink. The PV module is 
presented as a two diode model with the series’ and parallels’ resistances. 
The method is used to obtain the parameters of the module by using data 
provided by the manufacturer. To verify the accuracy of the proposed model, 
a comparison is done between the present results and those in a previous 
paper in which the single diode detailed model was studied. The same type of 
modules (PWX 500) is used to compare and validate results. The weather 
data as illumination and temperature are taken as input variables when 
current and Power are output variables. For the validation of the proposed 
model, it is used to power a resistive load by using a boost converter 
controlled by 4 Maximum Power Point Tracking method. The proposed 





Duty cycle method is very detailed and can be useful for PV researchers for their works 
Matlab/Simulink . : ; Me sot . 
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1. INTRODUCTION 

Electricity generated from fossil fuels which are based on oil and gas represents more than 60% of 
the world production. Until today, it remains the basis for electricity production. But it is also the first source 
of pollution. In addition to this, the traditional organization of the electrical system (production, transport, 
distribution) is seen over the decades, attenuated because of two key phenomena; population growth and 
change of the load quality which contains now more polluting devices such as computing devices 
for example. 

Renewable energy is needed nowadays as an alternative to fossil energy in depletion. Repetitive 
natural phenomena as sunlight, wind and wave offer to us a very good opportunity to recover the clean aspect 
of our planet. In recent decades, many researchers are working on the PV cell technology and its efficiency. 
The electric power generated by the PV module depends strongly on irradiance and temperature. 

The photovoltaic cell is the basic device which generates electricity from sunlight. Cells are grouped 
in series to increase voltage or in parallel to increase current. Modules and arrays are formed by the cited 
association of cells. The PV cell is presented by an equivalent circuit elaborated according to the 
fundamentals of semi conductors. The simplest model is an ideal equivalent circuit made of a current source 
in parallel with a diode as shown in Figure 1 [1] and which doesn’t take in consideration many parameters 
related to fabrication processes or to the materials 

A practical model includes a resistance Rs in series with the current source as in [2-6]. Rs is 
responsible for the reduction in output voltage [7]. However, this model exhibits serious deficiencies when 
subjected to temperature variations. A most practical model adds another resistance Ry in parallel with the 
diode as in [7-11]. Rp reduces the output current of the module because of the leakage current through the 
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parallel branch [7]. But the most accurate model is the Two Diode Rs Rp model [12, 13] as shown in 
Figure 2. 





Figure 1. PV cell equivalent circuits [1] Figure 2. Two Diode model of PV cell 


As manufacturers don’t give enough information on data sheets about the parameters which depend 
on weather conditions (illumination and temperature), some assumptions are necessary to establish a 
mathematical model of the PV cell and the PV module. After modeling the PV cell, the optimal power can be 
determined accurately with the mathematical model so that it matches with the experimental data. 
I(V) characteristic is a non-linear (Transcendental) equation with several parameters. Some of them are 
provided by constructors but have to be adjusted, and others must be computed in the present work as Rp, 
reverse saturation currents of the two diodes Is; and Is2, and their ideal factors A; and Az. The objective of our 
work is to determine the parameters and simulate the module with these settings in order to compare the 
power generated by simulation with the one obtained experimentally and delivered by the constructor. 
Then, the PV module is connected with a load via a DC/DC converter 


2. PRESENTATION AND MODELING OF PV MODULE 
2.1. One Diode PV model 

The red part of Figure 1 shows that the Kirchhoff law allows obtaining the output current I for the 
ideal model: 


I= lon Iq (1) 
Where In is the photocurrent, Ig the Diode current which is given by (2): 
oie 
Iq = Is [emma = 1| (2) 
Where V is the voltage imposed to the diode, Vr the thermal voltage and equal to (k.Tc/q), Tc is the cell 
temperature in Kelvin (K), k is the Boltzmann constant (1,381.10-°j/K), q is electron charge (1.602.10°°C), 


Ns is the number of cell connected in series, A the ideality factor and Is the reverse saturation current. 
The presence of the series resistance changes the diode current into the following form: 





(3) 


q(V+I.Rs) 
Iq = I, Jeanette = 1| 


The bleu part of Figure | shows the most practical One Diode model. The output current is also obtained by 
Kirchhoff law and written in (4): 


I =Ipn — 1a — Ip (4) 
Where Ip is the current of leak in parallel resistance Rp. 


According to the previous equation, the output current of the One Diode module with Ns cells in 
series is: 
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V+IRg 
Rp 





(5) 


a(V+I.Rs) 
T=Ip, —1, [ eanenre = 1| 


2.2. Two Diode PV model 
The most accurate model is presented in Figure 2. It is also called ‘’Two exponentials model’. 
It is given as: 








V+IRs 
Rp 


(6) 


q(V+L.Rs) q(V+I.Rs) 
_— Th — 154 | eA1Ns.k.T¢ -—1 — Is [ easnerte = 1] - 


3. DETERMINATION OF THE PARAMETERS 

Seven parameters must be determined Ipy, Is1,1s2, Rs, Rp, Ay, and Az. But as A, and A, are assumed, 
Is, and Is, are not determined and Iph is known, the parameters to be determined are only Rs and R, 
- The photocurrent under any conditions depends on irradiance and temperature and Iphyrer: 


lpn = = a Cones + kj. AT) (7) 


Where Ipnrer (A) is the photocurrent at STC and G is the Irradiance called also illumination (W/m’), 
Greg: Irradiance at STC = 1000W/m’, AT =T, —Terer (Kelvin), Toyer is the cell temperature at 
STC = 25+273=298K, k;is the coefficient temperature of short circuit current (A/K), provided by the 
manufacturer. 
- The reverse saturation current for each diode is given by the following expressions[19]: 








i —tecres) Cc (4 4") =) 

1s = Teper Oey (Soest) 7 (Ze) e xp ( Aqk YG *) (8) 
—tecrel) Tc (¢ . 

Ion = Iscref: exp (er) (Zz) e xp (ee 46) \(Z Toref -3)| (9) 


- Avis the ideality factor of the first diode according to the diffusion current densities. 
- Ag is theideality factor of the second diode according to recombination current density [17]. 
Most of papers on two diode model of PV cell assumed that A; = 1 and A2= 2 [7, 8, 13]. 

In the proposed model, A; = 1 and Az = 1.2 because when plotting P(V) characteristics, the best 
match is obtained with the practical model provided in manufacturer data sheet [14-18], and (6) is now 
rewritten at maximum power condition (according to data sheet): 

Pmax,ref (Vmaxref +Rs-Imaxref ) (10) 


= = 19, — Ta. - | 
ph,ref d1 d2 Rp 


vf = 
max,ref Vinanter 


Rg is varied to compute Rp in (6) which is implemented in Matalb Simulink environment. The 
iteration starts at Rs =0 and increases in order to force the computed max power point moving until it 


matches with the experimental max power point. The corresponding Rp is then deduced. Finally, only one 
pair (Rp, Rs) satisfies this condition and its expression at max power is written as follows: 


R= Vmax,ref +Rs.Imax,ref 
po qV+IRs) q(V+IRs) Pies (11) 
I —I eAiNs.k.T¢ —4|-] eAiNs.Kh.T¢ —4|——Seet es 
phref—ts1 2 Vinax.ref 








The PWX 500 PV module (49W) is used to validate the proposed model. The series resistance Rs is 
provided equal to 0.55Q. but Rp is not given as seen in Table 1. Results of the One Diode Model, previously 
given in detail in [19], are given now in table to be compared with results of the Two Diode Model and with 
the provided data at STC. 


Table 1. PWX 500 PV module (49W) characteristics at 25 °C, 1.54M, 1000W/m? 
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Provided parameters 


Adjusted parameters for 


Adjusted parameters for 





One Diode Model Two Diode Model 
P max 49 W 49 W 49W 
Imp 2.88 A 2.88 A 2.88 A 
Vinp 17V 17V 17V 
Isc 3.11A 3.11 A 311A 
Voc 21.8V 21.8V 21.8V 
Rs 0.55 QO 0.45 0 0.30 
Rp - 3100 200 Q 
Ai - 1.3 1 
Ao - - 1:2 
Isi - 4.116*10°A 1.782*10° A 
Is2 - - 9.075*10°A 
Noct 45 C° 45C° 45 C° 
Ki 1.3*10-3 (/K°) 1.3*10-3 (/K°) 1.3*10-3 (/K°) 
Ky -72.5*10-3(/K°) -72.5*10-3(/K°) -72.5*10-3(/K°) 
Ns 36 36 36 





4. SIMULATION OF THE PROPOSED MODEL 

Figure 3 represents the main system. Input variables are Temperature and Irradiance when output 
variables are Current delivered by the module and power generated to be used by load either in stand alone or 
grid connected mode. Figure 4 is a subsystem of the main system. It is an implementation of (6). It is used to 
simulate the proposed model by incrementing R, until matching A, with Ppp ex. The experimental data for 
maximum power at STC provided by the manufacturer of PWX 500 PV module (49W) are given in Table 1. 

To valid our model, the PV module is simulated under Matlab-Simulink. The general diagram is 
given in Figure 5. It comprises of a PV module and a resistive load in parallel. A diode is connected in series 
to prevent the reverse current flow. To make the voltage stable, a filter is connected before the load. 
A converter is used to boost voltage. Its gate is controlled by the output of MPPT (Maximum Power point 
Technique) bloc in order to get the max of the power delivered by the PV module. 

One of the most used methods of algorithms in MPPT is Perturb and Observe commonly called 
(P&O). It consists on making a perturbation in the voltage and observing the rate of change of power. 
It depends on the position of the optimum point to adjust the reference voltage [20]. A flowchart of P&O 
method is presented in Figure 6. 


oe 















Ea 


PV Module 


Figure 3. Presentation of the whole PV model 
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Figure 4. Detailed Two Diode model with Rp 
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Figure 5. General diagram of Photovoltaic system connected to a resistive load 
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Figure 6. The flowchart of P&O algorithm 


The duty cycle can control the DC/DC converter and maintain the voltage at the maximum. 
The diagram of Figure 7 taken from references [21] and [24] presents a part of our main system concerning 
the DC/DC converter built a high frequency IGBT switch. Authors give equations of inductor (L), capacitor 
(C) and the duty cycle (D): 


a Vint-D.T s _ Vin-D.T 2 _ Vout—Vin 
Ai,” AVo.R.f ? Vout 


L 


Where here, Ai, indicates inductor ripple current and AVg indicates capacitor voltage ripple, T is the time 
period (sec) and f the switching frequency (Hz) 











Figure 7. Circuit Diagram for Boost Converter [21, 24] 
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5. RESULTS AND ANALYSIS 

In a previous work [19], the One Diode Rp model was modeled and simulated under Matlab 
environment. The series resistance Rs was adjusted to 0.45 Q instead of 0.55 © provided by manufacturer 
and the parallel resistance was computed iteratively and equal to Rp = 310 Q (not provided). The module 
PWX 500 (49W) was used to validate results. In the present work, the same method is used to adjust Rs and 
compute Rp for the Two Diode Model. Adjusted and computed parameters are summarized in Table. 

Rs is iteratively increased and Rp is simultaneously computed. The iterative method gave Rs = 0.30 
and Rp = 200Q. These two values make the two diode model most representative for the chosen PV module. 
To simulate any other PV module, one can introduce in the (19) the respective experimental maximum power 
and then, the iterative method is used to determine the appropriate pair (Rs, Rp) which makes the model the 
most representative. 

After running simulation for some different values of Rs, I-V curves of are represented in Figure 8, 
and P-V in Figure 9. One can say in Figure 8, that the shape moves to the rectangular form when R, 
decreases. J,, and V,, are not affected by change of R,, so the fill factor which is given by (12) change only 
with Prp: 


FF =" (12) 


Voclsc 


This is in accordance with (11) where P,,, is with respect of Rs and Rp. 

When comparing the One Diode model curves and the Two Diode model ones as shown in 
Figures 10 and 11, one can say thatfor the first model, simulated curves match with experimental data 
at Rs = 0.45Q and Rp = 310Q. For the second model, simulated curve match with experimental data 
at Rs= 0.3Q and Rp = 200Q. 
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Figure 8. I (V) curves for different Rs Figure 9. P (V) curves for different Rs 
(Two Diode model) (Two Diode model) 























Figure 10. P (V) curves for different Rs Figure 11.1 (V) curves for different Rs 
(One Diode model) (One Diode model) 
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The proposed model can be used to simulate PV module at different levels of irradiance and 
temperature The I (V) and P (V) characteristics are respectively presented in Figure 12 and Figure 13, 
by varying irradiance from 200W/m? to 1000 W/m? at T =25°C. In Figure 13 and Figure 14, The I (V) and 
P (V) characteristics are respectively presented by varying temperature from 0°C to 60°C at G =.1000 W/m’. 
In Figure 14 and Figure 15, The I (V) and P (V) characteristics are respectively presented by varying 
temperature from 0°C to 60°C at G =.1000 W/m? 





Voltage (V) 


Figure 12. I (V) characteristics 
by varying irradiance at T= 25°C 





Power (WV) 


Voltage (V) 





Figure 13. P (V) characteristics 
by varying irradiance at T= 25°C 
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Figure 14. I (V) characteristics by varying 
temperature at G=1000W/m? 








Voltage (V) 





Figure 15. P (V) characteristics by varying 


temperature at G=1000W/m? 


The next step of the work is to use the proposed model in the topology presented in Figure 5. 
The one stage PV system adaptation is used to supply a resistive load. Since the one stage system is in 
continuous current mode, voltage and current at DC/DC converter (load side) are given by [22]. In the 
reference [23], the duty cycle D is taken as o and 77 is the efficiency of the boost converter: 


Vpv 
D 


1 


Ip=na- D Inv 


Vy, 





(13) 
(14) 


Figure 15 shows the simulation of power delivered to the load by varying irradiance at Is from 
400W/m? to 1000W/m? and from 1000W/m? to 400W/m? at 3s in the general diagram is given in Figure 5 
when using P&O algorithm. One can see that when varying irradiance, the power delivered to the load 
corresponds to the maximum power extracted to the module at each level according to the plots presented in 


Figure 13. 
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Figure 15 Power delivered to the load with P&O method 


6. CONCLUSION 

In this paper, an accurate model of PV cell is presented in details and simulated under 
MATLAB/SIMULINK environment. Elaboration of the presented model is in accordance with the 
fundamentals of semi conductors and the technology of the PV cell. The computed parameters of the PV cell 
are chosen because of their variation with irradiance and temperature. This model can be used as a tool to 
study all type of PV modules available in markets even if manufacturers don’t give all parameters. The series 
resistance Rs, is adjusted by increasing from zero and used to compute the parallel resistance in order to 
match the experimental maximum power point provided in data sheet with the computed one. In our case, 
a pair of (Rs =0.3Q , Re=200 Q) instead of(R;=0.55Q,and R, not provided) is obtained and then to simulate 
and validate the proposed model. In a previous work, a single diode Rp model was presented and this work 
about Two Diode Model is a new step. Then, the proposed model was simulated in a general diagram a 
resistive load via a continued bus. A DC/DC boost converter was used and controlled by an MPPT 
(Maximum Power Point Tracker) method. The chosen algorithm was P&O because it’s easily implemented. 
The accuracy of the proposed model is demonstrated by different plots which show that the power delivered 
to load match with the maximum power obtained from the PV module. In this work, the simplest form of 
load was used to explain and illustrate the proposed model. A possible future application of this model will 
be in grid connected PV system with a non linear load. 
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